Elliptic flow of thermal photons in relativistic nuclear collisions 
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We predict the transverse momentum (pt) dependence of elliptic flow of thermal photons for 
Au+Au collisions at the Relativistic Heavy Ion Collider. We model the system hydrodynamically, 
assuming formation of a thermalized quark-gluon plasma at some early time, followed by cooling 
through expansion, hadronization and decoupling. Photons are emitted throughout the expansion 
history. Contrary to hadron elliptic flow, which hydrodynamics predicts to increase monotonically 
with pt, the elliptic flow of thermal photons is predicted to first rise and then fall again as pt 
increases. Photon elliptic flow at high pr is shown to reflect the quark momentum anisotropy at 
early times when it is small, whereas at low px it is controlled by the much larger pion momentum 
anisotropy during the late hadronic emission stage. An interesting structure is predicted at interme- 
diate pt ~ 0.4 GeV/c where photon elliptic flow reflects the momenta and the (compared to pions) 
reduced v% of heavy vector mesons in the late hadronic phase. 

PACS numbers: 25.75.-q,12.38.Mh 



Experiments performed at the Relativistic Heavy-Ion 
Collider (RHIC) are providing evidence for the produc- 
tion of quark-gluon plasma (QGP) in nuclear collisions 
at ultra-relativistic energies. Key findings include strong 
anisotropic flow of all hadronic species 00,13 and a sup- 
pression of high-pT hadrons due to parton energy loss in 
the dense medium 0, jg". Signatures of direct photon 
emission [1, 0, H, , indicative of a hot early state, have 
also started emerging ^(J- The emphasis of the next 
generation of experiments will thus necessarily shift to 
a more precise determination of the properties of both 
the QGP and the subsequent hot hadronic matter. Hy- 
drodynamic flow, and in particular anisotropic flow in 
non-central collisions, provides strong evidence for the 
existence of a hot and dense initial state with thermal 
pressure Elliptic flow is generated very early, via 
the transformation of the initial spatial eccentricity of 
the nuclear overlap region into momentum anisotropics, 
through the action of azimuthally anisotropic pressure 
gradients. With the passage of time, the pressure gradi- 
ents equalize, and the growth of elliptic flow shuts itself 

off 13. 

Due to their strong interactions, hadrons decouple 
from the system late, typically when the temperature 
has dropped to values around 100 Me V 0>0|. In hydro- 
dynamic simulations, hadrons with large transverse mo- 
menta are thus emitted from those fluid elements which 
have the largest radial flow. This is supported by the 
hydrodynamically predicted p-p-dependence of the ellip- 
tic flow coefficient V2 and of the Hanbury-Brown Twiss 
radii Photons, on the other hand, are emitted at 

every stage of the collision, from the pre-equilibrium 
stage |H, the quark-gluon fluid, and the late hadronic 
matter. The thermal emission of pho tons from the QGP 
and hadronic phases is obtained |Fj| by integrating the 
thermal emission rate (which is strongly biased towards 



higher temperatures) over the space-time history of the 
system. As a result, high-px photons arise mostly from 
the hot early stage, where hydrodynamic flow is weak but 
the spatial eccentricity of the source is large, whereas 
hadrons are emitted when the temperature is low, the 
flow is strong and anisotropic, but the spatial eccentric- 
ity of the fireball has mostly disappeared. The elliptic 
flow of photons, especially at pr ^ 1 — 2 GeV/c, is there- 
fore expected to provide a glimpse of the early part of the 
expansion history when the fireball is in the QGP phase, 
complementary to the elliptic flow of hadrons. 

This last argument can be made independent of the va- 
lidity of hydrodynamics. For two of the leading photon 
production processes (quark-gluon Compton scattering 
and quark-antiquark annihilation for high energy pho- 
tons in the QGP, 7T7t — * pj for low energy photons in the 
hadronic phase) it is known that, for photon energies 
well above the rest masses of the emitting particle, the 
photon production cross section peaks very strongly for 
momenta close to that of the emitting particle jl4(. So, 
even without local thermal equilibrium, one expects the 
photon elliptic flow to track the momentum anisotropy 
of the photon-emitting particles at similar momenta. 

For Au+Au collisions at RHIC it is known that 
the hydrodynamic behaviour of elliptic flow begins to 
break down for mesons (baryons) above pr — 1.5GeV/c 
(2.3 GeV/c): instead of continuing to rise with pr, as pre- 
dicted by hydrodynamics |3J, the elliptic flow saturates. 
On the other hand, the systematics of hadron produc- 
tion with pt ^ 2GeV/c can be well described by quark 
coalescence jlfj. The observations then translate into a 
quark elliptic flow near hadronization that breaks away 
from hydrodynamics above pr — 0.75 GeV/c 0|- Thus, 
even though ideal fluid dynamics excellently describes 
the bulk of particle production at RHIC, viscous cor- 
rections become significant at transverse momenta above 
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1 GeV/c for quarks and gluons and above 2 — 3 GeV/c for 
hadrons. We expect these non-equilibrium features to be 
reflected by the photons, and the hydrodynamic predic- 
tion for photon elliptic flow presented here to be only an 
upper limit once pr exceeds 1 GeV/c. Nevertheless, the 
qualitative features pointed out below are generic and ex- 
pected to be robust against non-equilibrium corrections. 
The photon momentum spectrum can be written as 

EdN 1 /d 3 p = J [(...) exp(-p-u(x)/T(x))] d 4 x . (1) 

The square bracket indicates the thermal emis- 
sion rate from the QGP or hadronic matter. 
P p = (pt cosh Y, pt cos 4>, pt sin 4>, pr sinh Y) denotes 
the 4-momentum of the photon, while the 4- velocity of 
the flow field is w M = 7r(coshr?, v x (x, y), v y (x, y), sinh rj) 
(with 7t = (1— v^)' 1 / 2 , Vj,=v%+Vy), assuming boost- 
invariant longitudinal expansion. We use coordinates 
t, x, y, r), with volume element d i x = T dr dxdydr], 
where r = (i 2 — z 2 ) 1 / 2 is the longitudinal proper time and 
rj= taxih~ (z/t) is the space-time rapidity. The photon 
momentum is parametrized by its rapidity Y, transverse 
momentum px = (p 2 + P y v ) 1 ^ 2 , and azimuthal emission 
angle (f>. The photon energy in the local fluid rest frame, 
which enters as p-u/T into the Boltzmann and several 
other factors in the thermal emission rate, is given by 
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T 



[cosh(Y-77) - v T cos (4>-4> v )] , (2) 



where 4> v = tan~ 1 (v y /v x ) is the azimuthal angle of the 
tranverse flow vector. Eq. J5J shows that the azimuthal 
anisotropy (0-dependence) of the photon spectrum, con- 
ventionally characterized by its Fourier coefficients v n 
(where for equal nuclei only even n contribute at Y = 0) , 



dN(b) 

d 2 px dY 2npT dpr dY 
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[l + 2v 2 (p T ,b) cos(2</>) 



is controlled by an interplay between collective flow 
anisotropy and geometric deformation of the tempera- 
ture field T(x,y,r): It obviously vanishes in the ab- 
sence of radial flow, vt = 0, but for nonzero vt it can 
arise from an anistropic flow field or from an azimuthally 
deformed temperature field. Only if the transverse 
flow points radially (<j> v = 4> r = ta,n~ 1 (y/x)) and both the 
flow and temperature fields are azimuthally symmet- 
ric (vt{x, y) = vr(r) and T(x, y) =T(r) with r 2 = x 2 +y 2 ) 
can the ^-dependence in Eq. J2J) be made to disappear, 
by rotating in r = (x, y) — (r, 4> r ) such that 4> r = 4>. 

We emplo y t he boost invariant hydrodynamic code 
AZHYDRO jri| which has been used extensively to ex- 
plore hadron production at RHIC [l^ . We use stan- 
dard [lil Ht| | initial conditions for Au+Au collisions 
at — 200 A GeV, but extrapolated from the usual 
initial time To = 0.6fm/c to a 3 times smaller value 
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FIG. 1: Constant energy density contours for e — e q , eh, and 
et, along y(x=0) and x(y=0), for Au+Au at 6 = 7fm. 



of To = 0.2fm/c, assuming 1-dimensional boost-invariant 
expansion between these times. We do so in order to ac- 
count for at least a fraction of the pre-equilibrium pho- 
ton production at very early times |l8j : its contribution 
to the photon spectrum is important at large px, and 
will suppress v2 hoton there because very little transverse 
flow develops before 0.6fm/c. So our initial maximum 
entropy density in the center of the fireball for 6 = col- 
lisions is so = 3 x 1 17fm~ 3 = 351 fm~ 3 , corresponding to 
a peak initial temperature of Xo = 520MeV. The initial 
transverse entropy density profile is computed from the 
standard Glauber model parametrization 

We assume creation of a thermally and chemically 
equilibrated plasma at tq, and use the complete leading- 
order rate for QGP photon production from Ref. and 
the latest results for photon radiation by a hot hadron gas 
from Ref. 0- These rates provide good descriptions of 
si ngle photon data at SPS energies and recently at RHIC 
USES- We use tne equation of state EOS Q [13 which 
matches a non-interacting QGP to a chemically equili- 
brated hadron gas at T c = 164MeV, with energy densi- 
ties e q (T c ) = 1.6 GeV/fm 3 and e h (T c ) = 0.45 GeV/fm 3 in 
the subphases. Hadron freeze-out is assumed to happen 
at e { = 0.075 GeV/fm 3 [l^. 

Figure shows the changing spatial anisotropy of the 
fireball, for Au+Au collisions at a typical impact param- 
eter of b = 7 fm, by plotting cuts through constant energy 
density surfaces along the x and y-axes (i.e. in the reac- 
tion plane and perpendicular to it) for the three values 
e q , eh, and ef. After about 9fm/c the initial out-of-plane 
deformation changes sign (dashed and solid lines cross), 
as a result of faster expansion into the reaction plane. 
Figure shows the flow velocities along the same cuts, 
illustrating the developing anisotropy of the transverse 
flow field. Through most of the fireball interior the ve- 
locity along the x-axis is larger and rises more rapidly 
with increasing radius than along the y-axis. At hadronic 
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FIG. 2: Flow velocity along the constant energy density 
contours for e = e q (QGP phase, upper panel), eh (mixed 
phase, middle panel), and ef (hadronic phase, lower panel) 
for y(x=0) and x{y=0), for 200 A GeV Au+Au at b = 7fm. 



freeze-out, ef , the in-plane and out-of-plane velocity gra- 
dients have almost (but not fully!) equalized, and larger 
in-plane flow velocities exist mostly near the fireball sur- 
face. Figure El complements similar plots in Refs. 0,0] 
which show velocity profiles either at constant time or 
only along the freeze-out surface at 6 = 0. 

In Fig. [3] we show the thermal photon elliptic flow for 
Au+Au collisions at 6 = 7fm. We compare the overall 
V2(pt) with the individual elliptic flow coefficients asso- 
ciated with quark matter (^(QM)) and hadronic mat- 
ter («2(HM)) contributions to the thermal 7 spectrum. 
Comparison of ^(HM) with the elliptic flow of thermal 
7T and p mesons, ^(tt) and t>2(p), shows that at low 
(high) pt hadronic photon flow tracks the elliptic flow of 
ir (p) mesons. This can be understood from the single- 
photon spectrum (see, e.g., Fig. 17 in [HI): at low px 
hadronic photon emission is dominated by the processes 
p — ► 7T7T7 and 7T7t — > pr/ whereas at higher px collision- 
induced conversion of vector mesons into photons (such 
as np— > 7T7) takes over. In the present situation, with ra- 
dial flow effects on the spectra computed from hydrody- 
namics, this transition happens around pT~0.4GeV/c. 
At pT^0.4GeV/c hadronic photon elliptic flow thus 
tracks the vi of p and other vector mesons. Since the 
heavier vector mesons carry less elliptic flow than the 
lighter pions 0, this transition manifests itself in a re- 
duction of hadronic photon V2 around pT^0.4GeV/c 
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FIG. 3: (Color online) «2(pt) for thermal photons from off- 
central 200 A GeV Au+Au collisions at b — 7 fm. Quark and 
hadronic matter contributions are also shown separately, and 
the elliptic flow of ir and p mesons is shown for comparison. 

(see dashed red line in Fig. |3J). Due to the dominance 
of hadronic processes in the total photon yield in this px 
region, this structure survives in the total photon elliptic 
flow (solid red line). 

The elliptic flow of the QGP photons, t>2(QM), is small 
at low pr, peaks around pr~1.5GeV/c, and decreases 
again for large px- Its smallness at high pt and increase 
towards lower pt reflects the absence of transverse flow 
during the earliest, hottest stage of the QGP and its 
gradual buildup during the following cooler stages. That 
growth is cut off by the generic decrease of V2 as Br — > 
[22| . Its nonzero limit at px — (first observed in |23j for 
gluons) can be traced to the singularities of the Bose dis- 
tribution and emission rates at px — 0. Since the QGP 
photon rates are unreliable for pT<0.2GeV/c, the 
t>2(HM) curve in Fig. |3 should not be trusted in that do- 
main. Fortunately, there the total V2 (solid red line) is 
entirely dominated by the hadronic photon contribution. 

Although the hydrodynamically predicted elliptic flow 
of hadronic photons is almost everywhere much larger 
than that of the QGP photons, the hadronic contribu- 
tion to the photon spectrum is increasingly suppressed 
(by more than an order of magnitude below the QGP 
contribution) once px exceeds 1.5 — 2GeV/c. Hence the 
overall photonic V2, while larger than the pure QGP con- 
tribution, decreases for large px, too, approaching the V2 
of the QGP photons. At high px the photon elliptic flow 
thus opens a window onto the dynamics of the QGP, in 
spite of the larger elliptic flow of the hadronic photons. 

In Figure 01 we present the impact parameter depen- 
dence of the elliptic flow of thermal photons. The impact 
parameters are chosen to roughly correspond to collision 
centralities of 0-10%, 10-20%, 20-30%, 30-40%, 
40 — 50%, and 50 — 60% of the total inelastic nuclear 
cross section. As the impact parameter increases, the 
relative contribution of hadronic photons increases, too. 
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FIG. 4: Impact parameter dependence of the azimuthal 
anisotropy of thermal photons. 

Since its elliptic flow is larger than that of the QGP pho- 
tons, the 6-dependence of thermal photon V2 is predicted 
to be stronger than that of hadron V2- This can be used 
as an additional tool to isolate the QGP contribution 
to thermal photon emission, especially if methods can be 
developed |24| to subtract at high pt the hadronic contri- 
bution V2 (HM) from the thermal photon V2 , by exploiting 
its similarity to the elliptic flow of thermal p mesons. 

However, a quantitative interpretation of photon ellip- 
tic flow in the context of extracting dynamical informa- 
tion on the QGP evolution must also account for other 
contributions to high-p^ photon production. It has re- 
cently been argued that jets passing through the QGP 
contribute significantly to direct photon production at 
large pr and may even lead to negative photon ellip- 
tic flow in this region [25j]. It will be interesting to see 
how much of this effect survives once the so far neglected 
anisotropic collective flow of the medium is included. 

In conclusion, we have presented a first calculation 
of elliptic flow of thermal photons from ultrarelativis- 
tic heavy- ion collisions. The differential photon elliptic 
flow V2(pt) exhibits a rich structure, driven by the evo- 
lution of the system and the competing rates of photon 
emission from the quark and hadronic matter stages. Its 
impact parameter dependence differs from that of pion 
elliptic flow since the relative contributions of hadronic 
and quark matter emissions change with collision central- 
ity At high px, where ideal fluid dynamics is known to 
overpredict i>2, the elliptic flow of thermal photons from 
the hadronic matter tracks the flow of p mesons which 
can be measured independently. If this can be used to 
subtract the hadronic contribution to photon elliptic flow 
at high pt , a new window onto the early buildup of col- 
lective flow in the QGP will be opened. 
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